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ABSTRACT 

We report the results of an extensive VLT high-resolution imaging campaign of the 
central region of the giant globular cluster u> Centauri . More than 100,000 stars have 
been measured in the inner 9'x9' region of the cluster. On the basis of multiband 
Colour Magnitude Diagrams (CMD), we confirm the existence of multiple stellar pop- 
ulations along the Red Giant Branch (RGB). Moreover, thanks to the high-precision 
of this dataset, we conclude that the RGB does not present a smooth and continuous 
distribution, but shows a discrete structure: besides the metal-poor and the extreme 
metal-rich population (RGB-a), the existence of three metal intermediate populations 
is shown. 

Key words: techniques: photometric - stars: evolution - stars: Population II 
globular cluster: ui Cen. 



1 INTRODUCTION 

The understanding of the star formation history of the glob- 
ular cluster uj Centauri (NGC 5139) represents one of the 
most interesting challenges of stellar astrophysics, u) Cen 
■ is the most massive and luminous GC of the Milky Way 
$H 1 (M ~ 2.9 1O 6 M , Merrit et al. 1997) and it is surely the 
^ most peculiar in terms of structure, kinematics and stel- 
lar content. One of the most interesting peculiarities of ui 
Cen is its metallicity spread. The first sign of chemical in- 
homogeneity was found by Dickens & Woolley (1967) from 
the large colour width of the Red Giant Branch (RGB) in 
the colour- magnitude diagram (CMD). This result was later 
confirmed by a large number of spectroscopic works (Free- 
man & Rodgers 1975, Butler, Dickens & Epps 1978). 

In the last ten years, extensive spectroscopic surveys 
have been performed on large samples of giant stars (Nor- 
ris, Freeman & Mighell 1996, hereafter NFM96; Suntzeff & 
Kraft 1996, hereafter SK96), showing a multimodal heavy- 
element distribution. Recent wide- field photometric studies 
(Lee et al. 1999, Pancino et al. 2000 - hereafter POO, Hilker 
& Richtler 2000, Hughes & Wallerstein 2000, Rey et al. 
2004) have revealed the presence of an additional, well de- 
fined anomalous RGB (hereafter RGB-a). According to POO, 
this population contains approximately 5% of the cluster's 



stellar content and represents the extreme metal-rich end of 
the observed metallicity distribution. Moreover, the metal- 
rich giants appear to have different spatial distribution and 
dynamical behaviour with respect to the metal poor ones 
(Pancino et al. 2000, 2003, Jurcsik 1998, Norris et al. 1997). 

In this framework, we have started a long-term project 
devoted to the detailed study of the properties of the differ- 
ent stellar populations in this cluster (see Ferraro, Pancino & 
Bellazzini 2002) , based on a complete multi- wavelength pho- 
tometric survey and an extensive spectroscopic campaign for 
giant and subgiant stars, both in the optical and in the in- 
frared bands. Within this project, a number of results have 
been published on the photometric, chemical and kinematic 
properties of the RGB-a stars (see Pancino et al. 2000, 2002, 
2003; Ferraro et al. 2002; Sollima et al. 2004; Origlia et 
al. 2002, 2003) and the possible detection of its Sub Giant 
Branch (SGB, Ferraro et al. 2004). 

In this paper we present BVI photometry of more than 
100,000 stars in the central region of the cluster . In §2 
we describe the observations, the data reduction techniques 
and the photometric calibration. In §3 we present the CMD, 
describing the main characteristics of the different cluster 
populations. In §4 the photometric metallicity distribution 
for a large sample of red giants is presented and discussed. In 
§5 we show the detection of the RGB bumps of the different 
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populations of ui Cen. Finally, we summarize and discuss our 
results in §6. 



2 OBSERVATIONS AND DATA REDUCTION 

The photometric data were obtained with the FORS1 cam- 
era, mounted at the Unitl (UT1) of the Very Large Tele- 
scope (VLT) of the ESO (Cerro Paranal, Chile). All obser- 
vations were performed in Service Mode in 5 different nights 
on March, April and June 2002 (see Table 1), using the high 
resolution mode of FORS1. In this configuration the image 
scale is 0.1" pixel -1 and the camera has a global field of view 
of 3.4' x 3.4' . A mosaic of 9 (3 x 3) partially overlapping 
fields has been observed around the cluster center, covering 
a total area of ~ 9' x 9', (see Figure 0. 

After applying the standard bias and flat-field correc- 
tion, we used DAOPHOT II and the point-spread-function 
(PSF) fitting algorithm ALLSTAR (Stetson, 1987) to obtain 
instrumental magnitudes for all the stars detected in each 
frame. For each field, six different frames have been observed 
(4 long and 2 short-exposure) through the B, V and I filters. 
Repeated exposures have been combined in order to obtain 
a high signal-to-noise median image. The automatic detec- 
tion of sources was performed on the median frame adopting 
a 3(J treshold. The mask with the object positions was then 
used as input for the PSF-fitting procedure, that was per- 
formed independently on each image. As usual, the most 
isolated and bright stars in each field have been employed 
to build the PSF model (here a Moffat function of exponent 
2.5). For each passband, the magnitudes obtained from the 
long exposures were transformed to the same instrumental 
scale and averaged. The same procedure was followed for the 
short exposures, that were however used only to measure 
bright stars, which are saturated in the long exposures. The 
calibrating equations linking the instrumental magnitudes 
to the standard system were obtained from the comparison 
with the photometric catalog by Pancino et al. (2000, 2003). 
Finally, a catalog with more than 100,000 calibrated stars 
in the cluster has been produced. Our calibration has been 
compared with the photometric catalog by Rey et al. (2004) 
in the B and V passbands. The mean magnitude differences 
found are AB = 0.069 ± 0.066 and AV = 0.060 ± 0.054 , 
which are consistent with a small systematic offset in both 
passbands. 



2.1 Photometric errors 

The root mean square (rms) frame-to-frame scatter of in- 
strumental magnitudes computed for each star from re- 
peated exposures is a good indicator of the internal pho- 
tometric accuracy. Figure |2] reports the distribution of pho- 
tometric errors in different passbands as a function of mag- 
nitude. As expected, the errors significantly increase towards 
faint magnitudes, due to photon noise. 

We used this diagram in order to define a high-accuracy 
sample (hereafter HA sample) that contains only stars with 
the lowest rms. The behaviour of photometric errors has 
been modelled, as a function of magnitude, with analytical 
functions, adopting the curves that provide ay ms = 0.01 at 
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B=18, V=18 and 1=17. To construct the HA sample, only 
stars with errors below these curves have been considered. 
Figure[2]shows the selection performed on the global catalog: 
stars included in the HA sample are marked with heavy 
black dots. 
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Figure 1. Map of the region sampled by the observations. North 
is up, east on the right. Only stars with V < 16 have been plotted. 
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Figure 2. B (panel a), V (panel b) and I (panel c) photometric 
errors as a function of magnitude. Grey points represent measure- 
ments excluded from the high-accuracy sample. 



3 COLOR - MAGNITUDE DIAGRAMS 

Figures El and El show the (V, B — V) and (I,B ~ I) CMDs 
for the global sample of stars measured in uj Cen. Figure [5] 
shows the (V, B — V) CMD for the high-accuracy sample, 
where a more precise analysis of the sequences is possible. 
The main features of the presented CMDs are schematically 
listed below: 



(i) The CMDs sample the entire evolved population in the 
cluster, reaching the upper Main Sequence (MS) at V ~ 20; 

(ii) The anomalous RGB (RGB-a, as defined by POO), is 
also visible as a well separated population at the extreme 
red side of the main RGB structure; 

(iii) A complex structure of the intermediate RGB com- 
ponents (RGB-Mint (s)), which populate the CMD region 
between the RGB-MP and the RGB-a, can also be noticed; 

(iv) A narrow, well-defined Sub Giant Branch which 
merges into the MS of the dominant cluster population is 
visible in Figure at a magnitude level significantly fainter 
than the MS Turn-Off (this feature is described and dis- 
cussed in Ferraro et al. 2004). 

Apart from the anomalous SGB, the most surprising 
feature of the CMDs shown in Figures |3] 0] and |S] is the 
discrete structure of the intermediate RGBs. The complex 
morphology of the RGB has been already noted in several 
recent publications (e.g. Lee et al. 1999, POO and Rey et 
al. 2004). However, in all the above studies, the RGB-Mint 
population appears to uniformly populate the CMD region 
delineated by the dominant RGB-MP population and the 
recently discovered RGB-a. 

In the CMDs presented here instead a discrete, comb- 
like structure becomes evident, with a set of different in- 
termediate RGB-Mint (s) populating the region between the 
RGB-MP and the RGB-a. FigureEUshows a selection of RGB 
stars, done on the high-accuracy sample, at three different 
magnitude levels (V = 13, 13.5 and 14 respectively). The 
corresponding colour distributions are shown in three ad- 
jacent panels. As can be seen, five peaks can be identified: 
the two extreme peaks being the RGB-MP and the RGB-a, 
while three additional peaks can be distinguished between 
them (RGB-Mint (s)). 

This discrete structure of the RGB has strong impli- 
cations for the chemical evolution of u Cen, and suggests 
a star formation history characterized by bursts, as will be 
discussed in §7. 



4 METALLICITY DISTRIBUTION FUNCTION 

It is well known that the RGB colour of an old Simple Stel- 
lar Population (SSP) is mainly affected by the abundance of 
heavy elements and only to a lesser extent by the age of the 
population. The effect of age differences becomes smaller and 
smaller with increasing age, and for population II stars it can 
be neglected. Therefore, we used the RGB colour distribu- 
tion to derive a metallicity distribution function (MDF - see 
also Bellazzini et al, 2003). For this purpose, we selected a 
sub-sample of stars from the HA sample in the (V, B — V) 
plane, having V < 14.6, B -V > 0.75 (see Figure 0, in 
order to limit the analysis to the region in which colour is 
most sensitive to metallicity. We took special care in exclud- 
ing HB and AGB stars, which are easily identified from the 
high quality CMDs shown in Figures [3] [I] and [5] 

We then compared the selected RGB stars with a grid 
of ridge lines of galactic globular cluster (GGC), having 
different metallicities, taken from the sample of Ferraro et 
al. (1999 - hereafter F99). A metallicity estimate was de- 
rived for each star by interpolating the grid of templates in 
the (V, B — V) plane. A similar approach was adopted by 
Frinchaboy et al. (2002, hereafter Fr02) who compared the 
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Figure 3. (V, B — V) color-magnitude diagram for the global 
sample of ~90,000 stars measured in ui Ccn. 



Figure 5. (V, B — V) color-magnitude diagram for selcted stars 
in the high-accuracy sample (the HA sample). 




Figure 4. (I, B — I) color-magnitude diagram for the global sam- 
ple of ~80,000 stars measured in ui Cen. 



RGB star distribution from the Washington M, T2 photom- 
etry with theoretical isometallicity curves. Also Lee et al. 
(1999) and POO obtained empirical MDFs in co Cen by cal- 
ibrating, in terms of metallicity, the distance of each giant 
from the mean ridge line of the main metal-poor component 
(RGB-MP). Finally, Hilker & Richtler (2000) derived the 
photometric metallicity for a large sample of stars using vby 
Stromgren photometry and the related metallicity indices. 
In the following we analyze in detail the obtained MDF and 
compare our results with those of the authors quoted above 




Figure 6. RGB colour distribution at different magnitude levels 
(V = 13, 13.5 and 14 respectively), as marked on the left panel. 5 
different peaks are evident in the right panels, marked by arrows. 



and with the existing spectroscopic works of NFM96 and 
SK96. 



4.1 Basic assumptions 

In order to compare the RGB stars distribution with the 
mean ridge lines from F99, we need to assume a distance 
modulus and a reddening correction. In the following we 
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B-V 

Figure 7. The zoomed part of the high-accuracy (V, B — V) CMD 
in the RGB region. Stars selected to derive the MDF are marked 
with open circles. 

adopt the distance modulus by Thompson et al. (2001) 
(m — M)o = 13.65 ± 0.11 based on the observations of the 
detached eclipsing binary system OGLE GC17. Concern- 
ing the reddening and extinction coefficients, we assumed 
E{B - V) = 0.11 ± 0.01 (Lub, 2002), A B = LIE{B - V) 
and A v = 3.1E(B - V) (Savage & Mathis, 1979). 

In our analysis, we considered the global metallicity 
scale adopted by F99. Straniero & Chieffi (1991) and Salaris, 
ChiefR & Straniero (1993) showed that, when computing the 
isochrones of Population II stars, the contribution of the a- 
element enhancement can be taken into account by simply 
rescaling standard models to the global metallicity [M/H], 
according to the following relation 

[M/H] = [Fe/H] + log 10 (0.638 10 [Q/Fel +0.362) (1) 

Direct measurements of the cv-elements abundance in halo 
and disk field stars have shown a very well denned behavior 
as a function of [Fe/H], with a nearly constant overabun- 
dance ([a/Fe] ~ 0.3) for [Fe/H] < —1, and an evident trend 
with metallicity, which linearly decreases to [a/Fe] ~ 0.0 as 
metallicity increases (see Edvardsson et al. 1993, Nissen et 
al. 1994, Magain 1989, Zhao & Magain 1990, Gratton et al. 
1996). Since there is evidence (Pancino et al. 2002, Origlia 
et al. 2003) that the stellar populations in uo Cen have dif- 
ferent a-element enhancement levels, the global metallicity 
[M/H] is the most appropriate quantity to consider in the 
present case. 

4.2 Results 

Figure [H] shows the obtained MDF for 540 stars in the HA 
sample {panel a) and for 1,364 stars in the global sample 
{panel b), with V<14.6, B-V>0.75. The distributions are 
clearly asymmetric: beside the main peak, in both cases 



other four peaks are well visible. They can be associated 
with the following components: 

(i) The dominant metal poor population RGB-MP with 
[M/H] ~ -1.4; 

(ii) Three different metal intermediate components at 
[M/H] ~ -1.2 (RGB-Mintl), [M/H] ~ -0.9 (RGB-Mint2) 
and [M/H] ~ -0.7 (RGB-Mint3), respectively; 

(iii) The anomalous population RGB-a is clearly visible as 
a sharp peak at the extreme of the distribution (at [M/H] ~ 
-0.5). 

In order to check the validity of this approach, we com- 
pared our metallicity estimates with those obtained with 
the Ca Triplet survey by NFM96 (with the 41 stars in com- 
mon). To convert the metallicities presented by NFM96 into 
the global scale, we first converted the equivalent widths of 
NFM96 into the Rutledge et al. (1997) system, and derived 
the iron abundance in the Carretta & Gratton (1997) scale 
(see Section 5 in Rutledge et al. 1997). 

Then the derived metallicities have been converted into 
the global metallicity scale [M/H] using the eq. (1). Accord- 
ing to Ferraro et al. (1999), we assumed [a/-Fe] = +0.28 for 
stars with [Fe/H] C a < -0.8 and [a/Fe] = -0.35[Fe/H] C c 
for stars with [Fe/ H]cg > —0.8 (see F99). The comparison 
of NFM96 metallicities transformed adopting this procedure 
with our photometric metallicities is shown in Figure [§] As 
can be seen, the agreement is quite good. The solid line 
shown in Figure |U] represents the iso-metallicity relation. 
The bottom panel shows the residuals of the comparison. 
The nice agreement between the two metallicity determina- 
tions is confirmed within the errors, with a mean dispersion 
of (T[m/h] — 0.2 dex. 

The MDF shown in Figure |H] (a) is in good agree- 
ment with those presented by the other photometric sur- 
veys (Lee et al. 1999, POO and Fr02). The peak detected at 
[M/H] ~ —0.5 was not found by Fr02, probabily because, 
as pointed out by the author, their catalog is severely in- 
complete in the inner 6 arcmin, a region where the anoma- 
lous metal rich RGB-a seems to be concentrated (see POO 
and Pancino et al. 2003, hereafter P03). With respect to 
the spectroscopic works (NFM96 and SK96) and the photo- 
metric estimates by Hilker & Richtler (2000) we note that, 
while the global shape of the MDFs are similar, the metal 
rich components of the distribution presented in this pa- 
per are not reproduced. The lack of metal rich stars in the 
distributions of NFM96 and SK96 is probably due to two 
different effects: (i) the different mean radial distances of 
the samples: NFM96 and SK96 limit their analysis to outer 
regions of the cluster where the presence of RGB-a stars is 
less evident, (ii) the magnitude selection adopted by NFM96 
{V < 13), which tends to exclude RGB-a stars. Note that 
a few RGB-a stars are visible in the extreme red portion of 
the CMD selection boxes of SK96. Hilker & Richtler (2000) 
found a distribution similar to that of NFM96, with a contin- 
uous high-metallicity tail at [Fe/H] > —0.8. The different 
structure found at high-metallicities is probably due to the 
dependence of the Stromgren photometric metallicity index 
on CN abundance, which is not linearly dependent on metal- 
licity. 
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Figure 8. The derived MDF for the high-accuracy (a) and for 
the global sample (b). The significant peaks of the distribution 
are indicated in panel a and interpolated with a multi-burst en- 
richment model (panel b). 
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Figure 9. Comparison between the metallicities derived in this 
paper with those obtained by NFM96 for the 41 common stars 
(top panel). The residuals of the comparison are also shown (bot- 
tom panel). A solid line in both panels shows the iso-metallicity 
locus as a reference. 



4.3 Statistical significance of the peaks 

As quoted above, the MDF shown in Figure [S] suggests the 
existence of 5 peaks. Following the procedure described by 
NFM96, we assumed that the chemical evolution of a stellar 
system can be well represented by a simple model of chemi- 
cal enrichment (Searle & Sargent 1972, Hartwick 1976, Zinn 
1978, NFM96). This model considers the evolution of a cloud 
of gas having initial abundance Zq, which experiences is- 
tantaneous recycling of the ejecta from massive supernovae, 
and from which gas is expelled at a rate proportional to the 
star formation rate. Under these assumption the form of the 
abundance distribution is 



where 
fi(Z) 



h e 



f{Z)=Y,fi{Z) 



-{Z-Zoi)/(<Zi>-Zoi) 



< Zi > —Zen 



HZ) = o 



for Z > Z 0i 
elsewhere 



Where 1% represents the intensity of the i-th star forma- 
tion burst. 

In order to estimate the significance of the revealed 
peaks, we fit the derived MDF with this model, following 
the prescription of NFM96. Assuming [M/H] = log(Z/Z & ), 
we converted f(Z) into a function of the global metallicity 
f([M/H]) and smoothed it with a gaussian kernel of <r=0.05 
(equal to the width of the MDF bin, Zinn 1978). We there- 
fore performed a maximum likelyhood analysis of the data 
(see appendix B(c) of Morrison, Flynn & Freeman 1991). 
We treated the fit in terms of three free parameters for each 
burst: initial metallicity Zoi, average metallicity < Zi > and 
burst intensity U. 

The choice of the number N of bursts in the model must 
be done by balancing bias and variance of the fit: a small 
number of components leads to an uncertain fit, while large 
values of N lead to a more stable fit, that however might con- 
tain unnecessary components. For this purpose, we used the 
Bayesian Information Criterion (BIC - Schwarz, 1978). Let 
n (=32) be the bins of the MDF shown in Fig. 0(a), p (=3N) 
the free parameters of the fit and In the log- likelihood func- 
tion at the maximized parameters values 



log Ln — ^ ' log Pr 



3,N 



Where Prj t N represents the probability that the num- 
ber of stars whose metallicity lies in the j-th bin of the MDF 
be drawn from a N-bursts distribution. Then we choose N 
to maximize the quantity 



BIC(N) = I* 



- log n 



We observe that BIC(N) reaches its maximum value when 
N=4. This implies that the observed peaks in the MDF 
shown if Fig. |H| (a) are most probably reproduced by a four 
component model. 

The above analysis of the peak's significance is entirely 
based on a simple model of self enrichment. However, there 
is observational evidence suggesting that the case of to Cen 
can be much more complex than this (Norris et al. 1997). In 
particular, recent results cast some doubts on the hypothesis 
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Table 2. Relative frequency of the main populations of w Cen 





[M/H] 


% 


RGB-MP 


-1.4 


42 ±8 


RGB-Mint 1 


-1.2 


28 ±6 


RGB-MInt2 


-0.9 


17 ±5 


RGB-MInt3 


-0.7 


8±3 


RGB-a 


-0.5 


5±1 



that the RGB-a was the direct result of a self-enrichment 
process within u) Cen (Pancino et al. 2002, Ferraro et al. 
2002, 2004). In the following discussion, we will extend our 
analysis to all of the five components as they appear from 
the histograms in Figures and |H1 and we reserve to further 
address the statistical significance of each of them in the 
future. We will also treat and discuss the RGB-a separately, 
whenever needed. 

In order to estimate the relative frequency of each pop- 
ulation, we computed the area covered by each single burst 
distribution f(Zi). Since the selection performed on the HA 
sample could have corrupted the relative population ratios, 
we use the global sample. The results of the fit are plotted 
over the histograms in Figure |U(b) and summarized in Ta- 
ble where the peak metallicity and the relative frequency 
are listed for each considered population. We note that the 
extreme metal-rich population (RGB-a) represents ~5% of 
the stars of the whole system, fully in agreement with POO 
and P03. Regarding the metal poor population, its relative 
frequency is significantly smaller then the estimate of P03 
(~70%), while the Mint(s) populations represent a consis- 
tent fraction of the whole system with respect to what es- 
tabilished by P03 (~25%). This effect is most probably due 
to the contamination of the RGB-Mint 1 on the RGB-MP 
sample of P03, whose photometry did not allow a clear sep- 
aration between the two. 



4.4 Spatial distribution 

Using the RGB samples defined in the previous section, we 
investigated the spatial distributions of the sub-populations. 
In order to minimize the contamination between the metal- 
poor and Mints populations we limited our analysis to the 
brightest portion of the RGBs (V < 16, see Fig. EE). We 
computed the surface density distributions of each sam- 
ple using an adaptive kernel estimator algorithm (Fuku- 
naga 1972, Silvermann 1986, Beers 1991, Merrit & Trem- 
blay 1994). However, the smaller area covered by the pho- 
tometry presented here, and the subdivision in five com- 
ponents, drastically reduce the number size of each sample 
so that only a qualitative discussion can be done on the 
resulting spatial distributions. The comparison of the ob- 
served distributions reveals that both the Mints and RGB-a 
stars present a distribution that significantly differs from 
the RGB-MP sample's one. A two-dimensional generaliza- 
tion of the Kolmogorov-Smirnov test (Peacock 1983, Fasano 
& Franceschini 1987) gives a probability that the spatial 
distribution of the co-added (RGB-Mints + RGB-a) sam- 
ple and the RGB-MP one are drawn from the same parent 
distribution of less than 0.2 %. These results confirm those 
of POO, P03 and NFM96. Table El reports the location of 



Table 3. RGB centroids 
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201.703 


-47.4711 
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Mlntl 


201.694 


-47.4865 


-0.009 


-0.0154 


60 
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MInt2 


201.686 


-47.4850 


-0.017 


-0.0139 


64 


411 


Mlnt3 


201.744 


-47.4669 


0.041 


0.0042 


101 


244 


RGB-a 


201.691 


-47.4551 


0.012 


0.0160 


65 


133 
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Figure 10. Groups selection in the global (V, B-V) CMD. The 
different population samples are marked with different simbols. 
The continuous lines follow the distribution minima and are used 
to separate the groups. 

the centroids, the distance from the dominant population 
(RGB-MP) center and the number size of the sample for 
the different population of uj Cen. We observe that, while 
the centroid of the Mlntl sample is located close to the 
dominant MP component center, the three most metal rich 
populations (RGB-MInt2, RGB-MInt3 and RGB-a) seem to 
have different centroids. In particular, the RGB-a centroid 
appears displaced ~ 61" to the north of the dominant pop- 
ulation, in full agreement with P03. 



5 THE RED GIANT BRANCH BUMPS 

One of the most interesting features of the RGB is the so- 
called RGB bump (Iben, 1968). The location of the RGB 
bump in the CMD is mainly a function of age and metal 
content of a given stellar population. Because of the age- 
metallicity differences between the populations of ui Cen, 
one expects that this feature will appear in each RGB at a 
fainter magnitude as metallicity increases. The identification 
of the RGB bump is not an easy task because of the need 
of large observational samples (Crocker & Rood 1984, Fusi 
Pecci et al. 1990). 

To identify the RGB bump in our CMD we performed 
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the procedure described in the following. By computing the 
colour distribution of stars at different magnitude levels 
along the RGB in the global (V, B — V) CMD, we separated 
the five populations of ui Cen by following the position of 
the distribution minima (see Fig. EH. Then we located the 
RGB bumps by identifying the corresponding peaks in the 
differential RGB luminosity functions (LF). Figures HUl and 
I 111 shows the five RGB groups and their corresponding LFs, 
respectively. The bumps for the RGB-MP and RGB-Mint 
populations are clearly defined and marked with arrows, 
while the RGB bump of the extreme metal rich population 
(RGB-a) is much less evident and will be not considered 
in the following discussion. In Table [I] the observed Vbutup 
magnitudes for the four considered populations of ui Cen are 
listed. 

It is interesting to note that the location of the RGB 
bumps on our CMD does not appear as a "continuous and 
slanting feature", as it was described by Rey et al., (2004) 
but as a succession of distinct peaks. This is best put into 
evidence in Figure [T2l where the Hess diagram constructed 
on a zoomed area of the (V, B — V) CMD is presented. As can 
be seen, the bumps are clearly separated from each other, 
again pointing towards a set of well separated components 
in the RGB of w Cen. 

The absolute magnitudes of the RGB-bumps are plot- 
ted as a function of the photometric metallicities, obtained 
in Sect. 4, in panel a of Figure [*Hfl The result is compared 
to the empirical relation of F99 (obtained from the detec- 
tion of the RGB-bump in 54 GGC) and to the theoretical 
expectation by Straniero, Chieffi & Limongi (1997, hereafter 
SCL97), computed for Y = 0.23 and age= 16 Gyr. As can be 
seen (Figure^D, the location of the RGB bump for the four 
considered populations of u Cen is in agreement, within the 
uncertainties, with the predictions of both relations. This 
further supports the calibration of the derived MDF dis- 
cussed in Sect. 4 and the discrete nature of the RGB com- 
ponents in lu Cen. 

In panel (b) of Fig. 1131 the positions of the bumps are 
compared with isochrones in the plane M^" mp vs. [M/H] 
derived from the SCL97 models (Eq. 3 of F99). Neglecting 
the actual zero point of the resulting age-scale, that may 
strongly depend on the details of the theoretical models 
(see also Riello et al. 2003) , the distribution of the observed 
points over the isochrones grid suggests an age range of sev- 
eral Gyr between the MP and the most metal-rich Mint 
populations. While still within the observational uncertain- 
ties, the observed trend deserves some discussion, since it 
provides independent support to the results obtained from 
the analysis of Turn-Off stars with Stromgren photometry 
by Hughes et al. (2004), and with medium resolution spec- 
troscopy by Hilker et al. (2004). 

The derived age range (~ 6 Gyr) has to be considered 
as a strong upper limit because of (a) the unrealistically 
old zero-point of the age scale (18 Gyr) and (b) the ne- 
giection of the Helium enrichment which should be asso- 
ciated with the process that altered the mean metallicity 
from [M/H] ~ -1.4 to [M/H] ~ -0.7. The more metal rich 
bumps plotted in Fig. ll3l appear brighter than that expected 
if they were associated with population having the same age 
and He content of the MP population. Since both a younger 
age and an increase in He abundance (Y) should enhance 
the luminosity of the bumps, the observed trend may be in- 



Table 4. RGB bump parameters of lu Cen 



RGB bump 


Vbump 


[M/H] 


MP 


14.31 ± 0.03 


-1.4 ±0.2 


Mlntl 


14.42 ± 0.03 


-1.2±0.2 


MInt2 


14.57 ±0.03 


-0.9 ±0.2 


MInt3 


14.84 ± 0.04 


-0.7 ± 0.2 



terpreted either in terms of age and/or Y variations. In this 
context, several authors considered the effects of differences 
in helium abundance between the different populations of 
ui Cen (Ferraro et al. 2004, Bedin et al. 2004, D'Antona & 
Caloi 2004, Norris 2004). In particular, Norris (2004) sug- 
gested a helium enhancement of SY ~ 0.10 — 0.15 in order 
to reproduce the complex MS-TO morphology of lu Cen. 
SCL97 models predict 



Assuming that the MInt2 population corresponds to the 
"second population" of Norris (2004) (having A[M/H]=0.5 
with respect to the dominant MP population), the observed 
difference in the RGB bump level between these populations 
can be induced by an helium enhancement of AY ~ 0.09, 
once the effect due to the metallicity has been taken into 
account and assuming no age difference between these pop- 
ulations. This finding is not too far from what proposed by 
Norris (2004), but should be considered as an upper limit if 
we accept that MInt2 population is younger than MP one, 
as suggested by elementary concepts of chemical evolution. 
The most reasonable hypothesis is that both factors are con- 
tributing, according to the natural trend of chemical evolu- 
tion, e.g. successive generations of stars are younger than 
the MP population and more enriched in metals and He. 

In panel (c) of Fig. 1131 the differential Age- Metallicity 
Relation (AMR) obtained from panel (b) is displayed (A age 
= 0.0 is placed at age= 18 Gyr). This plot shows that the 
AMR derived from the RGB bumps is very similar to that 
obtained by Hughes et al. (2004) and it is fully compatible 
with the same "closed box" model adopted by these authors 
(see their Fig. 7). The above comparisons strongly indicate 
that there is full compatibility between the constraints ob- 
tained from the TO region of the CMD by Hughes et al. 
(2004) and Hilker et al. (2004). and the observed trend in 
the M^ mp vs. [M/H] described in the present paper. Both 
results suggest that the chemical enrichment of the u) Cen- 
tauri system in the range —1.4 ^ [M/H] —0.7 occurred 
on a timescale of a few Gyr (< 6 Gyr), the actual value 
depending on the unknown amount of Helium enrichment. 



6 DISCUSSION AND CONCLUSIONS 

We presented high precision BVI photometry for lu Cen, re- 
sulting from an accurate survey of the cluster's core. More 
than 100,000 stars have been measured, allowing a complete 
characterization of the RGB. In particular, the most ro- 
bust conclusion of our work is that the complex structure 
of lu Cen RGB does not show a smooth and continuous dis- 
tribution of populations, but it has a clearly discrete struc- 
ture. At least four, most probably five, separate components 
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Figure 11. Differential LFs for five RGB groups. The peaks cor- 
responding to the RGB bumps are marked with arrows. In the 
upper panel the groups selection in the (V, B — V) CMD is shown. 




B-V B-V 



Figure 12. Hess diagram for the RGB bump region marked in 
the left panel. In the right panel, the location of the RGB bump 
for the MP and the three Mints populations is marked with open 
circles. 



Figure 13. Panel (a): My of the RGB bumps as a function of 
global metallicity [M/H]. The solid line is the empirical relation by 
F99, the dashed line is the theoretical prediction by the Straniero, 
Chicfn & Limongi (1997) models at t=16 Gyr. Panel (b): Com- 
parison with the isochrones in the My vs [M/H] plane. The 
isochrones are derived from the same models that adopt a fixed 
He abundance of Y = 0.23 (see F99). Panel (c): Age-Metallicity 
Relation obtained from the bumps. The thick continuous line is 
the AMR predicted by the same "closed box" model adopted by 
Hughes et al. (2004), assuming constant star formation and yeld 
y = 0.005. 

have been identified. The discreteness and different nature 
of the components, clearly visible from the new CMD pre- 
sented (Fig. 0, has been further demonstrated by study- 
ing in detail: (i) the photometric metallicity distribution 
function, which shows five separate peaks (see Fig. |HJ and 
(ii) the identification of the RGB bump for each population 
(although some uncertainty remains on the RGB-a bump), 
which appear as distinct peaks if studied through a Hess di- 
agram (Fig. 1121 , and indicate an age difference of ~ 6 Gyr 
between the MP and the most metal-rich Mint population. 

Such a characteristic has one simple implication on the 
star formation history of u) Cen, if interpreted from the point 
of view of a simple self-enrichment scenario. It implies in fact 
that each population of uj Cen (except maybe the RGB-a, see 
below) can be associated with a different episode of star for- 
mation. Therefore, the chemical evolution of ui Cen should 
be viewed as a series of bursts, separated in time. However, 
a simple self-enrichment scenario may not be the best de- 
scription of the evolution of u Cen. We recall at this point 
the puzzling correlations between chemical and kinemati- 
cal properties found by Freeman (1985) and Norris et al. 
(1997). Also from the chemical point of view ui Cen shows 
some peculiarities which are not observed in normally self- 
enriching systems, i.e. dwarf galaxies and the field of the 
Milky Way. One such characteristic is the presence of CH- 
CN, Na-O, Mg-Al anomalies (Norris & Da Costa 1995a), 
which up to now has only been oserved in galactic glob- 
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ular clusters, u> Centauri and nowhere else. The second is 
the strong enhancement of s-process elements seen in both 
the metal-poor and intermediate components (Evans 1983, 
Norris & Da Costa 1995b, Smith et al. 2000, Vanture et al. 

2002) and the RGB-a (Pancino 2003) of w Centauri. This 
last peculiarity has only been observed in another stellar sys- 
tem, the Sagittarius dwarf galaxy, which is presently merg- 
ing with the Milky Way. Therefore, a simple self-enrichment 
cannot explain all the observational evidences gathered in 
the past, in particular when the RGB-a population is con- 
sidered. In fact it appears to stand out from the other popu- 
lations in terms of morphology, structure, chemical composi- 
tion and kinematics. The RGB-a main peak appears signifi- 
cantly dislocated from the main population's one. Moreover, 
the RGB-a chemical abundance, besides the higher iron con- 
tent, presents a smaller a-enhancement with respect to the 
other populations (Pancino et al. 2002, P03, Origlia et al. 

2003) while a large overabundance of heavy neutron capture 
elements remains high as for the other populations (Pancino 
2003). Finally, the mean proper motion of the RGB-a stars 
appears significantly different from that of the main clus- 
ter population (Ferraro et al. 2002). This last result have 
been questioned by Platais et al. (2003) who claimed that 
the different proper motion observed for the RGB-a is due 
to a spurious instrumental effect in the original proper mo- 
tion catalog by van Leuween et al. (2000). However, while 
the arguments claimed by Platais et al. (2003) were already 
considered and dismissed in the original publication by van 
Leuween et al. (2000) , they also have been strongly criticized 
by Hughes et al. (2004) who emphasised that no evidence 
of systematic trend has been detected in the proper motion 
residuals. Further support to the correctness of the proper 
motion measurements has been brought by Pancino (2003), 
who again did not found any significative residual trend as 
a function of either magnitude or color. 

Summarizing, the whole body of evidence points to- 
wards a complex scenario, where the different populations 
could have not only had a different chemical enrichment his- 
tory, but also a different dynamical evolution. It is clear that 
if one desires to accommodate all the observational facts in 
one self-consistent scenario, it is necessary to include not 
only a complex chemical evolution, but also a complex dy- 
namical history, possibily including minor mergers within a 
framework such as the "merger within a fragment" (Searle 
1977, Norris et al. 1997), or capture of smaller components. 
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